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1 Introduction 



In early eighties, Belavin and Drinfeld |BD| classified nonskewsymmetric clas- 



sical r-matrices for simple Lie algebras. It turned out that such r-matrices, 
up to isomorphism and twisting by elements from the exterior square of the 
Cartan subalgebra, are classified by rather unusual combinatorial objects which 
are now called Belavin-Drinfeld triples. By definition, a Belavin-Drinfeld triple 
for a simple Lie algebra g is a triple (Fi, T), where Fi, r2 are subsets of the 
Dynkin diagram F of g, and T : Fi ^ F2 is an isomorphism of Dynkin diagrams, 
which satisfies the nilpotency condition: if a G Fi then there exists k such that 
T'^~^{a) G Fi but T'^(a) ^ Fi. The r-matrix corresponding to such a triple 
is given by a certain explicit formula. This formula works not only for simple 
finite dimensional Lie algebras but in fact for any symmetrizable Kac-Moody 
algebra. 

In |s| , the second author generalized the work of Belavin and Drinfeld and 
classified classical nonskewsymmetric dynamical r-matrices for simple Lie alge- 
bras. It turns out that they have an even simpler classification: up to gauge 
transformations, they are classified by generalized Belavin Drinfeld triples, which 
are defined as the usual Belavin-Drinfeld triples but without any nilpotency 
condition. The dynamical r-matrix corresponding to such a triple is given by 
a certain explicit formula, which, as before, works not only for simple finite 
dimensional Lie algebras but in fact for any symmetrizable Kac-Moody algebra. 

This includes some well known examples: if T = id one gets Felder's dynam- 
ical r-matrix, and if g is of type An-i (i.e. the Dynkin diagram is an n-gon) and 
T is the rotation of the n-gon by an angle 27rfc/n where k is prime to n then 
one gets Belavin's elliptic r-matrix. 

G.Felder Q associated to every classical dynamical r-matrix, a remarkable 
system of differential equations called the Knizhnik-Zamolodchikov-Bernard (KZB) 
equations. For the Felder and the Belavin r-matrix these equations have a 
representation-theoretical interpretation. Namely, the KZB equations with the 
Felder r-matrix are satisfied by conformal blocks for the Wess-Zumino-Witten 
(WZW) model of conformal field theory on elliptic curves (or, equivalently, 
weighted traces of products of intertwining operators) and the KZB 

equations with the Belavin r-matrix are satisfied by conformal blocks for the 



WZW model twisted by the rotation of the Dynkin diagram p|, |KT]. It is 
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therefore natural to expect that a similar interpretation exists for all dynamical 
r-matrices. 

In this paper, we propose such an interpretation. Namely, for any Kac- 
Moody algebra and a (nondegenerate) generalized BD triple we show that 
weighted traces of products of intertwining operators, "twisted" by T, satisfy 
the KZB equations with the corresponding dynamical r-matrix from . 

We consider two cases: operators with values in representations from cat- 
egory O, and operators with values in finite dimensional representations for 
affine Lie algebras. In the first case we get the KZB equations for trigonometric 
dynamical r-matrices corresponding to generalized Belavin-Drinfeld triples for 
Kac-Moody algebras, and in the second case we get the KZB equations for ellip- 
tic dynamical r-matrices which are intermediate between Felder's and Belavin's 
elliptic r-matrices. 

In conclusion we would like to point out some directions of future research. 

First of all, it turns out that classical dynamical r-matrices mentioned above 
can be explicitly quantized. To obtain such a quantization has been an open 
problem, except for a few special cases, but we will present a complete solution 
of this problem in our next paper (joint with Travis Schedler). 

In another paper we plan to generali ze th e results of the present paper to 
quantum groups, following the ideas of |EV3], where it was done for the case 
Fi = r2 = r and T = Id. This will give quantum KZB equations, which are 
difference equations involving quantum dynamical R-matrices which are quan- 
tizations of the classical dynamical r-matrices that appeared in this paper. In 
the case of the Belavin R-matrix, these equations are the elliptic qKZ equations 
for the 8-vertex model, which play an important role in statistical mechanics. 



2 Classical dynamical r-matrices for Kac-Moody 
algebras 

In Q the second author associated a solution of the classical dynamical Yang- 
Baxter equation to every generalized Belavin-Drinfeld triple for simple Lie al- 
gebras. This construction easily extends to any Kac-Moody algebra. We recall 
this construction in this section. 



Preliminaries. Let A — [aij) be a symmetrizable generalized Cartan matrix 
of size n and rank I. Let (f), F, F) be a realization of i.e () is a complex vector 
space of dimension 2n — /, F = {ai, . . . a„} C f)* and F = {hi, . . . C f) are 
linearly independent sets and (a^, hi) — Uij. Let g = n_ © f) © n+ be the Kac- 
Moody algebra associated to A, i.e g is generated by elements ei, fi, i = 1, . . .n 
and [) with relations 

[cj, fj] = Sijhi, [(), [)] = 0, [h, a] = {ai,h)ei, [h, fi\ = -{ai, h)/^, 

together with the Serre relations (see [Q). Set b± = f) n±. Let t)' — 0Cft-i 
and let ( , ) be a nondegenerate invariant bilinear form on g. Recall that the 
restriction of ( , ) to () is nondegenerate (hence it defines a form on f)*, which 
we also denote by (,)), and that the kernel of the restriction of ( , ) to [)' is 
equal to the center c of g. Let A = A"'" U A~ C i)* be the root system of g and 
let Qa denote the weight subspace corresponding to a root a. For each a G A''" 
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fix bases \ . . . ei''"" and fa \ . . . fcf^"^ ^"^ of Qa and g-a respectively sucli 
that {ea\fa^) = Sij. The vector [ea\ fa^] is independent of the choice of i 
and we set ha = [ea\fa^] for any i = 1 . . . dim g^. Let p G f)* be an element 
satisfying (p, hi) = 1 for all i. Identifying t) with f)* via the form ( , ), we can 
regard p as an element of 1^ (p is well-defined up to adding a central element). 

Definition. A generalized Belavin-Drinfeld triple is a triple (Fi, T) where 
Ti, r2 C r and T : Fi ^ r2 is a inner product preserving bijection. 

Given a generalized Belavin-Drinfeld triple (ri,r2,T), we let F3 be the sub- 
set of Fi n F2 consisting of roots which return to their original position after 
applying T several times. It is clear that (Fi \ F3,F2 \ F3,T) is an ordinary 
Belavin-Drinfeld triple (i.e. T satisfies the nilpotency condition) and (F3, F3, T) 
is generalized Belavin-Drinfeld triple on which T is an automorphism of the 
Dynkin subdiagram F3. 

Set [ = {J2aeTi '^(^ ~ ^(^)))^ C t) and let l^o C f) be the orthogonal com- 
plement of [ in ^ with respect to the inner product on f). 

Definition. We say that the generalized Belavin-Drinfeld triple is nondegen- 
erate if the restriction of ( , ) to [ is nondegenerate. 
In the nondegenerate case, we have f) = I © f)o- 

Note that every generalized Belavin-Drinfeld triple on a Dynkin diagram of 

finite type or a connected Dynkin diagram of affine type is nondegenerate. In- 
deed, in a finite type case, [ is a real subalgebra of f), and the form on the real 
part of f) is positive definite, so [ is nondegenerate. A similar argument works in 
the affine case. Nondegeneracy is also the case when Fi = F2, because in this 
case T gives rise to a finite order (hence semisimple) orthogonal automorphism 
of f), and I = [)^. 

Let (Fi,F2,T) be a nondegenerate generalized Belavin-Drinfeld triple. It is 
convenient to choose an orthonormal basis {xj)j^i of f) with respect to ( , ) in 
such a way that 

[=0Cx„ 1)0 = ^Cxj, 
jell jeh 

for suitable disjoint subsets Ii,l2 C I such that / = /i U /2- 

Let i}i be the subspaces of t) spanned by h^.a G F^. 

The following Lemma is straightforward but important for the considerations 
below. 

Lemma 2.1. There exists a unique Lie algebra homomorphism B : n_ © f)i ^ 
n_®[)2 (resp. B'^ : n+(Si)2 ^ n+ffif}J such that B{fa) = /^(a), B(/i„) = 
ifae Fi, B{fa) =0 ifa^Ti (resp. B-^{ea) = eT-i(a), B'^iK) = /iT-i(a) 
ifa&V2, B-^{ea)=0ifa^T2). 

We will extend the homomorphisms B,B~^ of Lie algebras to the corre- 
sponding homomorphisms of their universal enveloping algebras. 
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Let (Ti), i G {1, 2, 3} be the set of roots a S A+ which are hnear combinations 
of simple roots from F^. Let be the subalgebra of g generated by 2a,Q-a, 
a € Ti. The map B restricts to an automotphism of gpa- For each a G (La), 
let Na be the order of the action of B on a, i.e B^°a — a but B'^a ^ a for 
1 < r < Na ■ Finally, it will be convenient to assume that fa ^ and ei* are eigen- 
vectors for B^", and set = ^- Note that B^°eL*^ = ei'^'^ei*^ 

Finally, let us define the Cayley transform of T on f)o- To do this, we need 
the following straightforward Lemma. 

Lemma 2.2. For any x d \)o, there exists a unique y G f)o such that for all 
a G Fi one has {a — Ta, y) — {a + Ta, x) . 

It is clear that y depends linearly on x. We will write y = Ctx. It is easy to 
check that the operator Ct : ()o ~> I)o is skewsymmetric. It is called the Cayley 
transform of T. 

Proposition 2.1 ([§). The function : I* ^ (g ® g)' 

oo 

is a solution of the classical dynamical Yang-Baxter equation 

i ^ 

+ [ri2(A),r"(A)] + [r^^ {X) , r^^X)] + [r^^ {X) , r^^ {X)] = 0. 

Remarks, i) In the expression for rxiX), the sum X^l^i ^"'^"'^''sa'' ^ ^'/q*"* 
is finite if a ^ (F3) and is an infinite series convergent to a rational function of 
e(">>) if a e (rg). 

Consider the special case when Fi = F2 = F, and T is an automorphism of 
the Dynkin diagram. Let N be the order of B (i.e. N is the smallest number 
divisible by all such that (6li*^)^/^° = 1). In this case, the formula for rxiX) 
can be written in the form 



2 ^ J 2 

^ „-i(a,A) 

Examples, i) When T = 1 one obtains Felder's trigonometric dynamical r- 
matrix 

r(A) = -^+ E ^cotanh(i(a,A))e(*)A/(*), 

a>0,i 

where fl = J2a ti^^^ ® f<^^ + Z"*^ ® "*) + X^i ® the Casimir element (see 



EV1|, Section 3). 
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ii) Consider g — sl{3) and let T be the automorphism exchanging the two 
simple roots. Then [ = Cp, so we can regard the element A S [* as a scalar. In 
this case, the dynamical r-matrix rxW is 

1 e"^-^ e"2A 

rr(A) = - - a;^ «) - ^ _ (eai /ai +6^2 8/02) ~ ^ , g-2A '^"i+"2 ® fa^+a^ 

~ "J ^_2A (•^"1 ® ^"1 •^"2 ® 6(12) ^ ^ p — 2A /"l+a2 ® 6qi+C(2 



Note that in this example the eigenvalue dai+a2 B is —1, which is the 
reason for the appearance of denominators 1 + e~^^. 

In the next section we give the representation-theoretic interpretation of the 
KZB equations associated to the dynamical r-matrix rT(A). 

3 Traces of intertwining operators 

3.1. Traces. For any f)-diagonalizable g-module V let V[X\ denote the sub- 
space of V of weight A G [)*. Let M\ be the Verma module of highest weight 
A and let vx S Ma [A] be a highest weight vector. Let Mj^ be the graded dual 
Verma module: = 0^A/[/i]* as a vector space and the Lie algebra q acts 

by 

x.a{u) — —a{x.u) V x £ Q, u E Ma, a G M^. 

Let t;^ G Afj^[— A] be the lowest weight vector satisfying {v'^,v\) = 1. Note that 
Ma and are irreducible for generic values of A. 

Recall the definition of the quadratic Casimir operator C: 

c=2p+j2x'^+2 y: Y.f^^- 

The operator C acts on Mx by multiplication by Aa := (A, A + 2p). 

Intertwining operators. Let A, /i G [)* and let y be a g-module from the 
category O. We will consider compositions of intertwining operators of the form 

$ : Ma ^ M^ ® F. 

The following lemma is well known (a proof can be found e.g in 
[|V|, |EFK1| ) : 



Lemma 3.1. Suppose that M* is irreducible. Then the map 

Homg(MA, M^ ® V[\ -p], $ {v;, $wa) 

is an isomorphism. 
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Given v G ^[A — /^], and when M* is irreducible, we will denote by <I>^ the 
unique intcrtwincr M\ ® V satisfying (w*, ^v\) — v. 

For any A £ f)* and any f)-semisimple g-module V we denote by A G End {V) 
the operator satisfying X\v[!y] = {^t'^)- 

Let G f)*. Consider the linear operator B : M^i defined by 

B(xVfj_i) — B{x)v^ for any x G Uri-. The following lemma is straightforward. 

Lemma 3.2. Let be such that {^\a) — {^,,Ta) for all a G Fi. Then for 
every x G C/(n_ © we have Bx — B(x)B and for every x G U{n+ © \)2) we 
have xB = BB'^{x). 

Let Vi, V2, . . . be g-modules from the category O, ui G Vi, . . . , G K- ho- 
mogeneous vectors of weights vi, . . .Vr respectively. Set v = ^Vi and consider 
(for generic /i) the composition 

■ ■ ■ ^'il'Be^ : M^, ^ M^_, ®Vi(^...(^Vr (3.1) 

where A G [*, n' G f)*, i^^',a) = {fi,Ta) for a G Fi (we regard [* as a subspace 
of f)* using the inner product). 



If 



(3.2) 



we can define 



pv 



-(A,A^)=Tr(CL,^_..._,_ 



.<S>l'-Be' 



B 



Vr. 



li V G I , the space of solutions of (3.2) in /i is an [* -principal homogeneous 
space. Note that it follows from the Kac-Kazhdan conditions (see p<K[) and 



from the fact that p G [ that for any fixed vi , . . . , Wr the composition ( |3.l| ) is 
defined for generic values of ^ in any [*-principal homogeneous space. In partic- 
ular, F'"^' " (A, /i) is a formal series in A whose coefficients are trigonometric 
functions of fi with values in the space {Vi® . . .® VrY . 



Set 

6b{\) = {TT\M_^{Be^)y' 
(a twisted version of the Weyl denominator). 

Lemma 3.3. We have 5b{\) = e^P'^^ UTse{r,)/B UtC^ " ^i'^e^^-f"^^)). 

Proof. Identify Un~ with M_p by u t-^ uv^p. Let C = Uq C Ui C ... 
denote the canonical filtration of Un^ . Since (L/„) is stable under Be^, we can 
replace Un^ by its graded when computing Tr|[/„- (Be^). By the PBW theorem, 

Gr C/n^ = Sn~ . By definition B acts nilpotently on fa '' unless a G {F3). Hence 
Tr|3n+(i?e^) = Trjgj^- (Be''') where rip^ = n~ fl qts- The Lemma now follows 
from 

aG(r3>/-B,t 
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and 

Tr 



S(C/<*>e...®Ci3«=-i/<")(-^^^) ^ _ g)(t)g-Ar„(aa) ' 



We put F"i-- '"'-(A,Ai) = (5i3(A)i^"i' -'"'-(A,^). This function is the main 
object of this paper. 

3.2. The KZB equations for traces. The following theorem is one of our 
main results. 

Theorem 3.1. The functionW^^'''' (X, fi) satisfies the following system of dif- 
ferential equations, for i — 1, . . . r; 

je/i ^ i>i 3<t ' (3.3) 

= ^ (A^-.,+i-...-.. - A^-..-...-.JF''^--'''-(A,m) 
These equations are called the KZB equations (see [^). 

Proof. For simplicity of notation we will write /Ui = /i — — . . . — i^i+i- We 
compute the function 

A{x, ^,) = Tr{^i\ . . . (q,,^^_^ - $;i;qM.^ . . . $;i^i?e^) 

in two different ways. On one hand, wc have 

A(A,^) = (A^,_, - A^Ji^"--'"-(A,/z). (3.4) 
On the other hand, using the relations 

(i^^e(^^)|,,^^_^ci>;^-ci>;;^(/K*)eW)|,^^^ 

-^M-Z'IM,, =-p\V.'^l\, 

we deduce that 

A(A,/i) ={ -25](/(*)e(*))|^_2p|^jF"--"'-(A,A.) 
+ Ai{X, ^i) + A2(A, Ai) + A3(A, ^) 

where 

Ai(A,^) = -2 ^x,|v.iv($;:; a:,|M,,_,$;:; ...<!';;'- se^), 

A3(A, A^) = -2 ^ il;l,lV(ci>- . . . ej)^^. . . . $-Be^). 
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Writing Ai as a sum of two equal terms: 

j&i 

and using the intertwining properties ^Xj = {xj^l + l®Xj)^ and (1 (8> ^xj) = 
$Xj — {xj 1)$ repeatedly in the first and second term of Ai{X, jj) respectively, 
we get 

jei \t<i t>i J 

- J2 (^(^i|M,, %\ ■ ■ ■ %^Be^) + Tvi^;\ . . . ^l^Xj\M,Be^)) 

(3.5) 

Now there arc two cases to consider, depending on whether j S 7i or j G l2- 
In the first case, by the cychcity of the trace, we have 

Tr(x,|M,,^;i • • • %^Be^) + TTi<Pl\ . . . $-^T,|M, Se^) 

=Tr{<i>;\...^;-{Bxj+xj)iM,Be^)- 

and 

Tt{^1\ . . . 1>;^{Bx, + x,)\mM) =2Tr($;; . . . ^I^x,\mM) 

d (3.6) 

=2-F--(A,,) 

where the differentiation is taken with respect to the parameter A. 

Let us now deal with the second case. It is easy to check using Lemma 3.2 
that for any a; e [-"-, one has 

Tr(a;|M,„ . . . %^Be^) + Tr($;:; . . . ^-;x\mM) 

= -Y^{CTx)iTv{^l\...^lrBe^). 

i 

Therefore, we obtain 

d 



^i(A,m) = ( - ^Xj\v,Xj\v, + ^Xj\v,Xj\v, - 2 E • 

^ t<i t>i jGh 

+ J2^mj2<^TXj\y,)F^^'-'^'^{X,f^) 

= ( - E^Ji^'^Ji^t + E^Ji^'^Ji^t ~ 2 E ^m-Q^ 



t<i t>i jeh 

+ E^^W E CtXj^v)f^^'--^^{X,h) 

(3.7) 

where in the last equality we used the skew-symmetry of "Ylij CrXj <Si Xj to get 
rid of terms with i = I. 
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We now compute A2{X, Using the intertwining property = (1 (81 
fa 1)^", we have 



f(t) 



= Tr(/i^i,^^ . . . + {fj^^^^ +... + f^^y^jF^^'- (A, 

(3.8) 



Hence by Lemma p.2| 



y CXj 00 \ 

^ i<i 1=0 j>i 1 = 1 ' 



(3.9) 



A similar computation (with moving to the right) shows that 



^ Ml Mi-l a 



(3.10) 

Adding (|j|) and ( plo| ) and using the relation 

'(5,Wy/JV„^^ if a G (Fg) and iV„|Z, 



else 



(3.11) 



we get 



A2(A, ^) + A3(A, ^) - 2 E(/^*^ei*^)v',^""-"'-(A, m) 

= -2fEE^iiv.(T^i("))iv,+EE^i?v,(^2("))i^.V"'---"^(^'/^) 



+ 2(EEC(?^3(a))|v,^----iA,M)+EEC(?^4(a))|v,)^----(A,A.) 

a,t j<.i a,t j'>i 

-2 E — 77T^T77-tW.^''"-''"-(A,A^)- 



(3.12) 



where 



(=1 
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Combining (^^, (|^ and ( ^.12D now gives the following relation for i^"!— '"-(A, ^): 
d \ 



where 

K{\) = y — 



A direct computation shows that 
_d_ 

' dx 



^ x~Tv\M^^{Be^) + {p + K{X))TT\M^^{Be^) = 0. (3.13) 



It is easy to deduce (3.3) from the above equations. 



3.3. The second order equation for traces. 

Theorem 3.2. The function F"i' " (A, fi) satisfies the following second order 
differential equation : 



q2 r \ 



F"i'-'"-(A,/i) (3.14) 



^ jeh 1 l,n=l 

where 



5t(A) ==^^^e-(^+^)("'^)(5VW ^5--ei*) +B-"e(,*) ®BV^*^) 

1 — Ct ^ 1 ~ Ct 



EL — 
- Xj (g) - Xi. 



jei2 

Proof. Consider 

A'ix,^l)^^T{^l\...^;^C\M,Be^)■ 

On one hand, we have 

and on the other hand, 

A\X, /.) =Tr ...$-( ^ x| + 2p + 2 ^ /l*^e(*0 |m, BeA 
^ jei a,t ' 

=AUA,A^)+A'2(A,A.)+A^(A,/i), 

where 

a; (A, y)^Y. • ■ • '^''^^mBe^). 

A',{X,^i) = 2T,{'^1\...'^I^p\mM) = 2^F----(A,m), 
A!,(A,A*) = 2^ Tr($- ...$-(/(*) e(*))|,,^Be^), 
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where the differentiation in the second equation is talccn with respect to A. To 
compute A[{X,fj,), note that, as in (|3.5|), ( |3.6D , 

= E ^^""•••'''-(^'^) + E Tr('i>;L. ■ --n^^M^Be')- (3-15) 



The second term on the r.h.s of (3.15) can be evaluated by the same method as 
in the derivation of the KZB equations: for all j G I2 we have 

Tr(<i>;i; x^%,^i?e^) =^^^x,|y,i^x,.|v.,F^-■-''■■(A,A.)^ (3.16) 

l,s 

We now compute A3 (A, fj,). We consider two cases 

Case 1: a ^ i^s)- By the intertwining property and the cyclicity of the trace 
again, we have 

Tv($;:;...<i>;:'-(/WeW)|M,se^) 

r 

= ^ UvM%\ ■ ■ ■ '^7^aUBe') + e-^"''^Tr($;:i^ . . . <i>--(e(*)i?(/(*)))|M.i?e^). 
1=1 

(3.17) 



Applying this equation repeatedly and using the relation (3.11), we obtain 

r OG 

(3.18) 



1 = 1 s=0 



Case 2: a e {T3). Applying ( |3.17| ) Na times yields 
Tr($;:i^...$;:'-(/(*)eW)|M.i?e') 

= E E e-^(-^)i?^(/.)|v,Tr(cI>-...$-e2)^^i?e^) 

(=1 s=0 

+ e"^°("'^)(?WTr(<i>;,; . . . 'f;^h^\M,Be^) 
+ e-^°(-^)0WTr(<i>;i; . . . $;:'-(/(*) eW)|M,Se^), 
from which it follows that 

TTi^l\...<fl''{fi''>e(^)^M,Be') 

= , (..) iE E e-^'-^^5^(/«)|v,Tr($-...$-ei*),,^i?e^) 

+ ei*)e-^°("-^)Tr($;:^^ . . . ^l^h^^M^Be')}. 

(3.19) 

The formula for A(j(A, a*) now follows from { ^.10\ ). Equations (|3.15D,(|3.16|), ( |3.19D 
imply the following second-order differential equation for i^''"!^- - (A, /i): 



2 ^ 9a;-^ dp 



Va^+ ^ 5t(A)^,«^„)f"--'"-(A,m) 



Ln=l 



(3.20) 
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where 



In particular, 



1 d_ 

l_6((*)e~JVc(«,A) dha 



d 



H{\) )Tr|M_„(i?e^) - iApTr|M_„(i?e^) 



2 dx) ' dp 



which, together with (|3.13D and (|3.20|) yields (|3.14D. 



3.4. Diagonalization of the KZB and the second order operators. 

Denote by Kj{X) the differential operators appearing on the left hand side of 
the KZB equations, and by D(X) the second order operator appearing in The- 



orem 3.2. These are operators on the space of functions of A with values in 
(Vi ® ... ®Vrf ■ It is known Q that Kj commute with each other. Besides, it 
can be shown that the operators commute with D (in fact, this is also clear 
from the discussion below). This gives rise to the problem of simultaneous diag- 
onalization of these operators. More precisely, the problem can be formulated 
as follows. 

Fix a weight G l^. Fix a generic point ^ G 1*, and consider the space of 
formal series 

W,.^^ e(^^«+'')(Fi ® ... ® V;)M[[e-(""^)]]. 

It is clear that the operators Kj , D act naturally in this space and are up- 
per triangular with respect to the natural ordering. The problem is to find a 
(topological) basis of W^,^ in which these operators are diagonal. 
The following proposition provides such a basis. 

Let Bi be homogeneous bases of Vi., and B(y) be the set of collections 
(ill, Wr) of vectors Vi G Bi such that the sum of their weights is v. 

Proposition 3.1. For generic the functions F^^'--'"''(A, ^^^"^ v + ^^—Y^Uiai), 
where rii > 0, and (ui, ...,Wr) run through B, form a common topological eigen- 
basis of the operators Kj, D in the space W^^^. 

This proposition follows immediately from the theorems of this section: the 
fact that the listed functions form a basis is obvious, so the only thing to be 
shown is that they are eigenfunctions, which was shown above. 



3.5. Quantum integrable systems associated to generalized Belavin- 
Drinfeld triples for simple Lie algebras. In the case when the Lie algebra 
g is finite dimensional, one can define other differential operators which commute 
with Kj,D. 

Namely, if Z is any element of the center of U{q) then there exists a unique 
differential operator Dz on 1* with values in End((Vi ® ... ®) Vr)[i'^]) such that 

Dz¥ ^ SB{X)TTi<i>;\ . . .<fl^Z\M,Be^)- 

For example, Dc = D — {p, p). 

It is easy to see that Dz^Zi — Dz^Dz^^, so the operators Dz form a com- 
mutative algebra. It is clear that these operators also diagonalize in the basis 
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of the previous section. Thus, we get a "quantum integrable system", whose 
eigenstates are the functions F"'^'--'"''. 

In the special case Ti = T, T — id, this system is a generahzed trigonometric 
Calogero-Moser system considered in [EFK2]. 



4 Classical dynamical r-matrices with spectral 
parameter 

Applying the construction of Section 2 to an (untwisted) affine Lie algebra 
g and using the evaluation map ev^ : — > 0, one can obtain solutions of the 
classical dynamical Yang-Baxter equation with spectral parameter. This is done 
as follows. 



Let g be a simple complex Lie algebra, and let Q ^ Q[t, t^^] © Cc © Cd be the 
associated affine Kac-Moody algebra, where c is the central element and d is 
the grading element. The commutation relations in g are: 

[c,g]=0, [xt",yf'^]^[x,y]r+"'+nSn,-^c, [d,xt"]=nxt^ Va;,yeg. 

Recall that the Cartan subalgebra of g is ^ = () Cc © Cd. Then I)* — 
\)* © CAo © C6 where Aq, S are defined by (Aq, i)) = (Ao, d) = {d,t)) = {S, c) = 
and ((5,9) = (Ao,c) = 1. Under the standard bilinear form on f), c and d are 
orthogonal to t) and we have (c, d) = 1, (c, c) = (9, d) = 0. 

The root system of g is A = (A + Zd) IJZ*5. The root subspace corresponding 
to a + kS is spanned by eat'' if a G A+, f-at'' if a G — A+ and equals f)t*'' if 
a = 0. The system of positive roots is A+ = A+ U U (A + N5). 



We will consider a twisted version of the affine Lie algebra g. Let g be the 
dual Coxeter number of g and set e — e~ . Consider the automorphism 7 = 
Ad(e2"''/f) of g. We have7(e„) = el"lea,7(/a) = and7||, = Id. Let g^ 

be the subalgebra of g consisting of all elements a{t) + Xc + ^d satisfying a{et) = 
7(a(t)). The elements eatl"l+™9, a;,i"f,c,9 for a e A+,m G Z 

and (xi) an orthonormal basis of [}, form a C-basis of g-y. The proof of the 
following lemma is straightforward. 

Lemma 4.1. The map 4> '. Q-y ^ Q defined by 

./.(e^il^l+^s) = e^i", 0(/at-l"l+™^) = fat"', Hd) =9d + p 

0(x,i'"'') - (m^O), ^{x^) = x, - {p, x,)- , 0(c) = - 

9 9 

is a Lie algebra isomorphism. 

Let (ri,r2,r) be a generalized Belavin-Drinfeld triple for g^. Let I be the 
subalgebra of i) of elements x such that {a,x) = {T{a),x) for a £ Fi. Let 
I — I ni). It is clear that [ contains c and 9, so I = Cc © Cd © [. 
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For any z G C* let ev^ : Q[t,t~^] ^ g be the evaluation map defined by 
evz(a;i") — z"x for all a; G g, n £ Z. 

Fix a complex number r with positive imaginary part. Let A = A + 2iTiT5/ g. 
Define 

rriX, z) = (ev^ ® evi)(0"^ ® (/)"^)(rT(A)), 

where z e C*, A G 1*. 

Remark. Although the evaluation maps are not defined on 5, this definition 
makes sense, since d occurs in in a combination c®d + d®c, and evz (c) = 
for any z. 

It is clear that Tt is a Laurent series in z whose coefficients are meromorphic 
functions on 1*. 

Proposition 4.1. The series rT{X, z) : [* — > is convergent in a nonempty 

annulus, and extends to a meromorphic function on [* x C*. Moreover, this 
meromorphic function satisfies the dynamical Yang-Baxter equation with spec- 
tral parameter: 

J2x[''±r^^iX, - ) A,i3(A, fl) + .f ) A,i2(A, fl) 

^ ' ^X^ 23 dXi Z3' dXi Z2 

+ [r'^iX, '-^)yHX, ^)] + [ri3(A, r'^X, ^)] + ^ (A, r'^X, ^)] = 0, 



where Xi is an orthonormal basis of [) . 

This proposition follows easily from Proposition 2.1. 

Now we would like to compute friX, z) explicitly. For the sake of simplicity 
we will restrict ourselves to triples of the form (Fi = r2 = T,T) (i.e T is 
an automorphism of the Dynkin diagram): these triples give rise to elliptic 
dynamical r-matrices. 

Remark. More general triples give rise to partially trigonometric and par- 
tially elliptic r-matrices. 

Let T be an automorphism of the Dynkin diagram F of g of order N . As 
in Section 2, let B be the lift of T to an automorphism of g of order N . By 



Lemma 4.1, B defines an automorphism of g-y. Note that the action of B on 
preserves the principal gradation, i.e the exists a unique automorphism (3 

of g such that B{xt™) = (3{x)t"^ for any xt"^ G g^. Furthermore, c and d are 

_B-invariant, and [) is i?-stable. Like before, we will choose an orthonormal basis 

Xi of f) which is compatible with I = f)-^ and [-'-. 
Let 

00 

j=-oo 

be the Jacobi theta function. For brevity we will not write the dependence on 
T explicitly. Introduce the functions 

e{w - u)e'{Q) e'{u) 



e{w)e{u) ' ' 9{u) ' 



where 9' is the derivative of 0. 
Let z = e^'^Ws. 
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Proposition 4.2. The function rxiX, z) : I* x C* ^ (g ® g)' is equal to 

a>0 1=0 
a>0 1=0 

('/lere a runs through positive roots of g). 

The proof of this proposition is by a direct calculation. 

Thus we see that the above construction gives a dynamical r-matrix with the 
number of dynamical parameters equal to d — 1, where d is the number of orbits 
of T on the Dynkin diagram. 



Examples, i) When T = 1, this yields (up to a gauge transformation) Felder's 
elliptic r-matrix, as shown in [EV1|, Section 4.6. 

ii) Let g = sl(n) and let T be the rotation of the Dynkin diagram by 2-Kk/n. 
Note that in this case (and in this case only) d = \ and we obtain a non- 
dynamical r-matrix with spectral parameter. It is easy to check that it is equal 
to Belavin's classical elliptic r-matrix with modulus r (see ]BD| |). 



5 Twisted traces of intertwiners for affine Lie 
algebras 

In this section, we apply the same procedure as in Section 3 in the case of an 
affine Lie algebra g, but we consider finite-dimensional (evaluation) modules Vi 
rather than modules from the category O. We keep the notations of Section 4. 

The Lie algebra has a triangular decomposition = f) © n~ where n+ 
(resp. n-) is the subalgebra spanned by eatl"l+'"», ra > 0, , x^t^f, 

m > (resp. spanned by to < 0, eail"l+'"f, Xii'"^, m < 0). 

This decomposition allows us to define highest weight Verma modules and dual 
Verma modules in the usual way. We put fl = 11+ |Ao — 2(k+g) ^ ^"^"^ denote by 
M^^k the Verma module M^. 



Finally, we define evaluation representations of g^. Let be a finite-dimensional 
highest weight g-module with highest weight vector wq. For any A £ C let 
z^^V[z, z^^\ denote the g-evaluation module, with g-action given by 

a;r.(wP(z)) = xvz'Tiz), c.{vP{z)) = 0, d.{vP{z)) = vz-^P{z) 
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for allx e Q, V e V, P{z) G z-^<C[z, z-^]. Let z-^V-y[z, z-^] C z-^V[z,z-'^] 
be the subspace of all (Laurent) polynomials v{z) satisfying v{ez) — e^'^'" p I ^ v{z) . 
Then z~^V^ [z, z^^\ is a g^-module (which is isomorphic to the usual g-evaluation 
module). From now on, we simply write z^'^V^[z^"'^] for z~'^V^[z^"'^]. 

Proposition 3.1 admits an analogue in this situation (see e.g [ EFK1 |,[^]) : 

Proposition 5.1. Let /i G ()* and fc G C and let V be a finite- dimensional 
Q-module. Suppose that M* ^ is irreducible. Then for each v G V[X — /i] there 
exists a unique g-y intertwiner 

l>;:,fc(z) : Mx^k ^ M^,fe®z-^y[z±l] 
where A — g ''^'2{k+g)^'^ such that 

Here ® denotes the completed tensor product. 

Let Vi, . . .Vr be finite-dimensional g-modules, and let ui G Vi, . . . G Vr be 
homogeneous vectors of weight vi, . . .Vr respectively. Set v = '^Vi, and assume 
G . We will consider composition of intertwining operators : 

^l-,,-...-u,.,ki^i) ■ . ■ $;:fe(z.)Se^ : AV,fc ^ M^,_,,k^z-^'Vi[zf'](g,. . .^z-^^Vr[2 

where A G (f)*)^, M € i)*, where {il',a) = {fi,Ta) for a G Fi, and A; = 
_^ 2{v^,^i-v^^i-...-Vr)-{v^,v^) ^ rp-j^jg composition lives in a certain completion (see 



[|EFK1[ , §3) of the space 

Hom(AV^fe, M^_„^fe ® Fi ® . . . (g) K-) ® z'^' . . . z-'^-C[zf\ ... , z^^]. 

Let us fix some t G C. For AGf)*setA = A + 2TTiTS/g. Consider the trace 
function 

Fr'-"'-(A,/i, z) = Tr ($- • ■ • *;:fe(^0i3e^) 

where z = {zi,. . . , z,.). Finally, let 

Sb{X) ^ (Tr^MojBe^))'' 

and set 

F^--'''-(A,M,z) =fe(A)Ffe(A,M,z). 
Our main results in the affine case are 

Theorem 5.1. The function F]^^'"'""^ (A, /i, z) satisfies the following system of 
differential equations for i — I, . . .r : 

"■^rr^^^T^^fc'^'"""^^'^'^) = Il2;,iy,^F;;;-"'-(A,/i,z) 

g uzi uxj 
3>i j<i 

(5.1) 
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These equations are called the KZB equations (see [^). 



Theorem 5.2. The function F"^' " '"'■(A, /i, z) satisfies the following second or- 
der differential equation : 



(5.2) 



where 

St{X, z) := {ev, ® evi){(f>-^ ® 0-1)(^t(A)), 
and S't(A) is the function defined in Theorem 3.2 (for q). 

The proofs of these theorems are completely parallel to the proofs of the 
theorems of Section 3. 

Remark 1. If Fi = r2 = F then ^^(A, z) can be expressed in terms of 
elliptic functions. 

Remark 2. In section 3 we saw that the differential operator 
D = Ai — S't(A) (where A; is the Laplacian of 1) on functions with values 
in the space of weight zero under [ in some representation of the Lie algebra 
can be included in a quantum integrable system. A similar statement holds for 
the operator Dr = A[ — ^^(A, 1), whose coefficients are elliptic (if F^ = F). 
To obtain other operators commuting with Dt, it is necessary to apply the 
construction of Section 3.5 to central elements of a completion of U{q) at the 



critical level k = —g, as explained in [ EFK2 |. In the cas e T — 1, this gives a 



generalization of the elliptic Calogero-Moser system (see [ EFK2 |). 

Remark 3. Although we considered only untwisted affine algebras, the 
results can be easily generalized to the twisted case. 
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